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ABSTRACT 

We present the results of rest-frame ultraviolet spectroscopic observations of a sample of 14 z ~ 3 
star-forming galaxies in the SSA22a field. These spectra are characterized by unprecedented depth in 
the Lyman-continuum region. For the first time, we have detected escaping ionizing radiation from 
individual galaxies at high redshift, with two of the 14 objects showing significant emission below 
the Lyman limit. We also measured the ratio of emergent flux density at 1500 A to that in the 
Lyman-continuum region, for the individual detections (C49 and D3) and the sample average. If 
a correction for the average IGM opacity is applied to the spectra of the objects C49 and D3, we 
find /i5oo//900,corr,C49 = 4.5 and fi.»snn/ km.mrr m = 2.9. These n umbers imply Lyman-continuum 
escape fractions at least as large as that presented in lSteidel et all |2001), measured from a compos- 
ite spectrum of 29 Lyman Break Galaxies (LBGs). The average emergent flux-density ratio in our 
sample is (fi .'snn/ fann,rnrr )= 22, implying an escape fraction ~ 4.5 times lower than inferred from the 
ISteidel et al \ composite spectrum. If this new estimate is representative of LBGs, their contribution to 
the metagalactic ionizing radiation field is J^(900) ^ 2.6 x 10"^^ erg s"^cm^'^IIz~^sr~^, comparable 
to the contribution of optically-selected quasars at the same redshift. The sum of the contributions 
from galaxies and quasars is consistent with recent estimates of the level of the ionizing background 
at z ^ 3, inferred from the H I Lya forest optical depth. There is significant variance among the 
emergent far-UV spectra in our sample, yet the factors controlling the detection or non-detection of 
Lyman-continuum emission from galaxies are not well-determined. Specifically, the two objects with 
detections are not differentiated in a consistent manner from the remainder of the sample in terms 
of their spectroscopic or photometric properties at wavelengths longer than the Lyman limit. There 
are also differences in the average emergent far-UV spectra of current and previous samples used to 
estimate the escape fraction from star-forming galaxies. Because we do not yet understand the source 
of this variance, significantly larger samples will be required to obtain robust constraints on the galaxy 
contribution to the ionizing background at z ~ 3 and beyond. 

Subject headings: galaxies: high-redshift — intergalactic medium — cosmology: observations — diffuse 
radiation 



1. INTRODUCTION 

The metagalactic ionizing background determines the 
evolution of the physical state of the intergalactic 
medium (IGM). The ionization rates of H I and He II 
are controlled by both the amplitude of the radiation 
field and its spectral shape between 1 and 4 Rydbergs. 

^ Based, in part, on data obtained at the W.M. Keck Observa- 
tory, which is operated as a scientific partnership among the Cal- 
ifornia Institute of Technology, the University of California, and 
NASA, and was made possible by the generous financial support 
of the W.M. Keck Foundation. 



These properties of the ionizing background reflect the 
total and relative contributions of quasars and star- 
forming galaxies. It is possible to estimate the na- 
ture of the background from the quasa, r proximity ef- 
fect l|Baitlik et al.ll988tlSco'tt et alJ200(il) . and also from 
the evolution of optical depths of H I and He II in the 
Lya forests of both those species, which are probed 
by the spec tra of bright quasars (ShuU et al. 200^ 
Bolton et alJ l2005l 12001 iMcDonald fc Miralda-Escudl 
200lUFan et al.ll2005(l . Such studies indicate that, while 
the comoving ionizing emissivity decreases between z — 3 
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and z — 5, the decrease in the comoving space density of 
quasars is much steeper, and that, by z 5, quasars 
cannot dominate the ionizing backgrou nd llFan et al.l 
|E)01, 2005; McDonald & Miralda-Escudl l2001() . Even 
at z ~ 3, the mean H I optical depth, and the mean 
and variation in the ratio of He II and H I column den- 
sities, 7V(He II)/A^(H I), may indicat e a substantial con- 
tribution from sta r-forming galaxies ijBolton et al.ll2005l 
|T006; Shu U et al.l l20041. Since Lyman-continuum radia- 
tion from galaxies appears to constitute an increasingly 
important ingredient of the ionizing background at red- 
shifts beyond z ^ 3, and may domina te the radiation 
field at the epoch of H I reionization (! Fan et all 1200 It 
ISpergel et al.ll200'3j) . estimating the galaxy contribution 
directly represents a fundamental goal for observational 
cosmology. 

One inherent complexity in such an estimate stems 
from the unknown value of fesc, the escape fraction 
of ionizing photons from the interstellar medium of a 
star-forming galaxy. Knowledge of this escape frac- 
tion is required to convert the star-formation rate den- 
sity of a population of galaxies into a comoving ioniz- 
ing emi ssivity. a crucial component of models of reion - 
ization l)Furlanetto et al.ll2004t iLoeb fc Barkanall200H) . 
Theoretical models of the propagation of ionizing radi- 
ation through the ISM of disk galaxies similar to the 
Milky Way pr edict an escape fraction of fesc < 0.1 
ijPove fc ShullllT99l iDove et alJl2000|) . though this frac- 
tion may increase for high-redshift galaxies with active 
star-formation rates, in which supernova explosions in- 
crease the ISM porosity ( Clarke & Oey 2002) and large- 
scale outflows may clear out chi mneys through whic h ion- 
izing radiation can easily escape (jFuiita et al.l2003l) . The 
model predictions for fesc depend on the detailed struc- 
ture and dynamics of the ISM, driven by the effects of su- 
pernova explosions and stellar winds. Currently, we have 
only fairly crude data on the prope rties of the ISM in dis- 
tant galaxies ijPettini et al.ll200H IShaplev et al.ll2003j) . 
with minimal spatially-resolved information. 

Another approach is to measure directly the emergent 
Lyman-continuum radiation from galaxies. Hopkins Ul- 
traviolet Telescope and Far- Ultraviolet Spectroscopic Ex- 
plorer (FU SE) observations of a small sample of nearby 
starburst gala xies mostly yielded upper limits for ff,,,r of 
less than 0.10 ('Leithere r et alJll995l:lHurw itz et"a nil997t 
iDeharveng et al. 2001). Only recently, Bcrgval l et al.l 
^2005) reported the first FUSE detection of Lyman- 
continuum emission from a local dwarf starburst galaxy, 
Haro 11, which implies fesc ^ 0.04 — 0.10, and will pro- 
vide a useful low-redshift baseline for comparing with the 
observations at high redshift. 

There has been some degree of controversy over the 
value of fesc in galaxies at z ~ 3. Esca ping Lyman- 
contin uum flux was apparently detected bv lSteidel et al.l 
l)200lD in a composite spectrum of 29 Lyman Break 
Galaxies (LBGs) drawn from the high-redshift tail of the 
z 3 LBG selection function. If the ratio of Lyman- 
continuum to 1500 A flux density observed in this com- 
posite spectrum is typical for the entire LBG sample, it 
implies fesc ~ 0.1, after correction for the average IGM 
opacity and the estimate that onl y 15 — 20% of 1500 A 
photons escape from typical LBGs l)Adelberger fc Steidell 
l200Qf) . However, the spectrum may not be representa- 
tive of typical LBGs, as it reflects the mean emergent 



spectrum for LBGs with bluer than average rest-frame 
UV colors and, therefore, objects less affected by dust 
extinction. Other attempts to measure escaping ioniz- 
ing radiation f rom LBGs have yielded only upper limits. 
Giallongo ct al. ( 2002) targeted two objects drawn from 
the bright end of the LBG luminosity function, obtaining 
an upper limit on the e scape fraction t hat is 4 times lower 
than the detection in iSteidel et all lj200H) . This limit 
may not be as stringent if the two individual sightlines 
pass through higher than average H I optical depth in the 
Lyman-continuum region - i.e. application of the average 
IGM correction factor for z ~ 3 may not necessarily be 
correct in individual cases. iFernandez-Soto et al.l p003) 
fit the broad-band spectral energy distributions (SEDs) 
of 27 galaxies at 1.9 < z < 3.5, including the effects 
of intrinsic Lyman-limit absorption, and intervening Ly- 
man series and Lyman limit absorption, and found on 
average fesc < 0.04. It should be noted that there are 
significant uncertainties associated with this method, be- 
cause of the use of the broad HST/WFPC2 F300W filter 
to probe the Lyman limit at z ~ 2 — 3. For galaxies to- 
wards the low redshift end of the sample, galaxy flux 
from above the Lyman limit falls in this filter, while for 
galaxies towards the high redshift end of the sample, the 
level of flux within the filter will be modulated primarily 
by the effects of IGM absorption, providing only crude 
information on the flux level just below the Lyman limit. 
Using a narrow-band filter more finely tuned to probe 
just below the Ly man limit for two galaxies at z ^ 3, 
ITuoue et a1J ^Q0) also found only upper limits to the 
escape fraction, but imaging with greater depth will be 
required to t est if the limits are in contradiction to the 
ISteidel et all pOOL) detection. 

Because of the increasing IGM optical depth at higher 
redshift, it is not possible to measure escaping ionizing 
radiation directly from galaxies much beyond z 3. Un- 
derstanding how fesc relates to other galaxy properties 
that can be measured at z '-^ 6 is required for assessing 
the contribution of galaxies to the ionizing background 
at the epoch of reionization. Therefore it is crucial to re- 
solve the controversy over the escape fraction from galax- 
ies at z ~ 3. In an effort to address this question, we 
have undertaken a program of deep spectroscopy of a 
sample of 14 LBGs, drawn from the heart of the LBG 
redshift selection function ((z) = 3.06). The mean and 
range of UV colors in this sample are representative of 
the total sample of LBGs, in con t rast t o the bluer than 
average sample of ISteidel et all l)2001l) . Spectra were 
obtained with the upgraded LRIS spectrograph on the 
Keck I telescope, which provides unparalleled sensitiv- 
ity in the 3500 — 5000 A wavelength range ( Steidel et a^ 
120041 . The depth in the Lyman-continuum region for 
individual spectra in our sam ple is comparable to that 
in the composite spectrum of ISteidel et al.l l|2001l) . In 
order to correlate the physical properties of the ISM in 
these galaxies with the degree of Lyman-continuum es- 
cape, we have also obtained sensitive observations red- 
wards of Lya for our sample, covering strong low- and 
high-ionization interstellar absorption features at twice 
the resolution of typical LBG discovery spectra. 

For the first time, we present the detection of escap- 
ing ionizing radiation for individual galaxies at z ~ 3. 
Considering our two detections and 12 non-detections, 
we find significant variation in the escape fraction from 
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galaxy to galaxy, and an average escape fra ction which is 
signifi cantly lower than the one presented in lSteidel et al.1 
()2001l) . Samples an order of magnitude larger are re- 
quired to understand the cause of variation among the 
emergent far-UV spectra of LBGs in terms of stellar pop- 
ulations, morphology, orientation, and ISM physical con- 
ditions, and to generalize to the global contribution of 
galaxies to the ionizing background. In § 2, we present 
our observations and reductions, while § 3 contains the 
main empirical results. In § 4, we discuss the emergent 
far-UV spectra of galaxies in our sample, and conclude in 
§ 5 with the implications for the contribution of galaxies 
to the ionizing background. A cosmology with f2„j — 0.3, 
JIa = 0.7, and h = 0.7 is assumed throughout. 

2. OBSERVATIONS AND REDUCTIONS 

2.1. Dte'p Spectroscopic Observations 

We targeted galaxies in the SSA22a field at 
RA=22:17:34.2, DEC=-h00:15:01 (J2000) for deep spec- 
troscopic observations. This field contains 146 photomet- 
ric LBG candidates, 57 of which were spectroscopically 
confirmed at the tim e the deep observations were planned 
ijSteidel et alJl200^ . Our primary targets for deep ob- 
servations were 14 bright galaxies with 23.0 <TZ< 24.5, 
twelve of which had redshifts previously determined from 
their original discovery spectra. With the fairly high S /N 
observations described below, we also measured redshifts 
for the remaining two galaxies. The redshift distribution 
is characterized by an average of (z) — 3.06 ± 0.12, and 
ranges from z = 2.756 to z = 3.292. Nine of the galax- 
ies in our sample have redshifts that place them within 
the significant overde nsity contained in t he SSA22a field 
at (z) = 3.09 ± 0.03 ijSteidel et al.llT99^ . However, the 
SSA22a field was simply chosen in order to maximize the 
number of bright {TZ < 24.5) targets with known spec- 
troscopic redshifts that could all be placed on the same 
multi-object slit mask, without regard for the presence 
of the redshift spike. 

The data were obtained during five separate observ- 
ing runs in the interval 2000 October - 2002 August, 
using the Low Resolution Imaging Spectrometer (LRIS; 
lOke et allllMl on K eck I in its newly- commissioned 
double-beamed mode ijSteidel et alJl200^ . Conditions 
were photometric, with seeing ranging from O'.'S— 0'.'9 for 
the majority of the observations, and never worse than 
I'.'O. As there is not currently an atmospheric dispersion 
corrector for LRIS, all but a few of the exposures were 
taken at elevations within 30° of zenith to minimize the 
effects of differential refraction in the blue. The total 
range in airmass spanned from 1.06 to 1.27. Two differ- 
ent slit mask configurations were used in the course of 
data collection, containing slits with 1'.'2 width. The slit 
lengths varied from 10" to ~ 50", with a median of 
~ 20" . The masks had different sky position angles but 
significant overlap in targeted objects, with 9 objects in 
common between the setups. On the 2000 October ob- 
serving run, the slit mask had a sky position angle of 
9 = 172°. The "d560" dichroic beam splitter was used 
to direct light bluer than ~ 5600 A towards the LRIS 
blue channel, while longer wavelength light was sent to 
the red channel. The blue light was dispersed by a 400 
lines mm^^ grism blazed at 3400 A, leading to a disper- 
sion of 1.70 A/pixel on the engineering-grade Tektronix 
2K X 2K UV/AR-coated CCD, while the red light was 



dispersed by a 600 lines mm^^ grating blazed at 7500 A, 
yielding a dispersion of 1.28 A/pixel on the red side Tek- 
tronix 2K X 2K CCD. Most slitlets had complete wave- 
length coverage spanning from the atmospheric cutoff at 
3200 A to redder than 7500 A. On all subsequent ob- 
serving runs (2001 June, 2001 July, 2002 June, and 2002 
August) a slit mask with a sky position angle of = —50° 
was used. The "d500" dichroic beam splitter was used to 
split the incoming light beam at 5000 A. The blue light 
was dispersed by the same 400 lines mm~^ grism blazed 
at 3400 A, while the red light was dispersed by a 600 lines 
mm~^ grating blazed at 5000 A. For the final two observ- 
ing runs (June 2002 and August 2002), the new science- 
grade Marconi 4K x 4K CCD mosaic was used for the 
blue-side observations. The new blue camera has smaller 
pixels, yielding a dispersion of 1.07 A/pixel when used 
together with the 400-line grism. For this second setup, 
most slitlets had complete wavelength coverage spanning 
from the atmospheric cutoff to redder than 7000 A. 

Our typical observing strategy consisted of acquiring 
a series of 1800 second exposures simultaneously on the 
red and blue sides. The telescope was dithered slightly 
(~ 1 arc second) between exposures in order to sample 
different parts of the detectors. We obtained 6 hours of 
total integration time for both red and blue side obser- 
vations in 2000 October, ~ 3 hours in 2001 June, ^ 5 
hours in 2001 July, ^ 2 hours in 2002 June, and ~ 8 
hours in 2002 August. Not all of the data acquired were 
included in the final stacked spectra. For red side spec- 
tra, data from 2001 June, 2001 July, 2002 June, and 2002 
August were included; for blue side spectra, data from 
2000 October, 2001 June, 2001 July, and 2002 August 
were included. At wavelengths below 4000 A, only 2002 
August blue side data were included in the final spectra, 
as the sensitivity of the science-grade blue detector in 
this wavelength range is a factor of 2 better on aver- 
age than that of the engineering-grade detector.^ The 
resulting total integration is ~ 8 hours below 4000 A, 

22 hours between 4000 and 5000 A, and ^ 17 hours 
above 5000 A. For a minority of the sample, only a subset 
of the this integration time was available. Tabled lists 
the TZ magnitudes, G — TZ and J7„ — G colors, Lya emis- 
sion and interstellar absorption redshifts, and total LRIS 
integration times in the wavelength ranges A < 4000 A, 
4000 < A < 5000 A (blue side), and A > 5000 A (red 
side), for the 14 objects targeted with deep spectroscopy. 

Spectroscopic flat fields for the red side data were ob- 
tained at the end of an observing sequence using an in- 
ternal halogen lamp. While flat fields were not used 
to correct blue side data taken with the engineering- 
grade detector, twilight flat flelds were obtained during 
the August 2002 observing run using the science-grade 
blue detector. The science-grade detector is character- 
ized by a wavelength-dependent flatfield pattern that 
becomes especially pronounced below 4000 A. Fluctu- 
ations resulting from imperfect correction of the fiat- 
field pattern constitute one of the largest sources of sys- 
tematic uncertainty in our measurements of flux in the 

^ Red side data from 2000 October were not included because 
of only partial overlap in wavelength coverage with subsequent ob- 
servations, and because of problems with background subtraction 
caused by internal reflections. Blue side data from 2002 June were 
not included because the necessary twilight flat fleld exposures were 
not obtained. 
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Lyman-continuum region, as we attempt to quantify in 
section 13.11 An observation of internal arc lamps (Hg, 
Ne, Ar, Kr, Xe for the October 2000 run, and Hg, Ne, 
Ar, Zn, Cd for all subsequent runs) was also obtained 
at the end of each observing sequence for the purpose 
of wavelength calibration. In order to flux-calibrate the 
deep mask spectra, a spectrophotometric standard star 
was observed at the end of each night through a long 
slit of width I'.'O, and with the same dichroic, grism, and 
grating combination as the deep mask. 

2.2. Data Reduction 

The data were reduced using IRAF tasks, with scripts 
designed for cutting up the multi-object slit mask im- 
ages into individual slitlets, flatfielding, rejecting cos- 
mic rays, subtracting the background, averaging individ- 
ual exposures into final summed two-dimensional spec- 
tra, extracting to one dimension, wavelength and flux- 
calibrating, and shifting into the v acuum frame. Thes e 
procedures are described in detail in lSteidel et al.ll)2003D . 
There was one substantive difference in the data reduc- 
tion procedures applied to the deep spectroscopic obser- 
vations, relative to those used for typical LBG discovery 
spectra. Rather than basic redshift identification, the 
goals of the deep observations include sensitive absolute 
measurements of flux in the Lyman-continuum region 
and accurate measurements of the depths of interstellar 
absorption lines at longer wavelengths. To achieve these 
goals, we took considerable care with the process of back- 
ground subtraction, which affects both of these types 
of measurements. Specifically, during the background- 
subtraction procedure, we made sure to exclude from the 
fit to the background all pixels within ± I'.'O of the ob- 
ject continuum location at each dispersion point. This 
technique eliminated most of the oversubtraction that 
is observed in LBG spectra with bright continua. Also, 
while most slits were significantly longer than 15" , we did 
not use more than this length to estimate the sky back- 
ground level in the Lyman-continuum region. A limit 
of 15" provided ample slit real-estate for estimating the 
background level, while preventing artifacts at large slit- 
distances from biasing the fit near our targets. 

The above procedures were applied to data from in- 
dividual observing runs. After the data from each run 
were fully processed to one-dimensional, wavelength- 
calibrated, flux-calibrated, vacuum-frame spectra, we 
combined the spectra from different observing runs for 
individual objects. For a given object, the spectra ob- 
tained on different observing runs ranged in S/N by as 
much as a factor of ~ 3 — 4, with the 2002 August spec- 
tra providing the highest S/N. Therefore, in construct- 
ing the final summed spectrum for each object at wave- 
lengths longer than 4000 A, we weighted the spectrum 
from each observing run by the square of its S/N. As de- 
scribed earlier, below 4000 A, only data acquired in 2002 
August with the science-grade detector were included in 
the final spectra. The final S /N per resolution element in 
the continuum redwards of Lya ranges from 5 — 25 for 
our sample, with a median value of ^ 12. With 1'.'2 slits, 
the resolution was determined by the seeing and result- 
ing object sizes, which we measured along the slit spatial 
dimension to be smaller than the slit-width in almost all 
cases. In the final combined spectra, we estimate typical 
resolution elements of 4 A and ~ 7 A on the red and 



blue sides, respectively. 

For 9 of the 14 galaxies, both Lya emission and average 
interstellar absorption redshifts were measured, whereas 
only interstellar absorption redshifts were measured for 
the 5 galaxies with no detectable Lya emission. Three 
galaxies were characterized by double-peaked Lya emis- 
sion. The average offset between Lya emission and in- 
terstellar absorption redshifts for this sample is consis- 
tent w ith what is observed in a much larger sample of 
LBGs (iShaplev et a l. 2003). While the current sample of 
spectra represent a significant improvement over typical 
LBG discovery spectra, in terms of signal-to-noise, they 
are still not deep enough for the robust measurement 
of the wavelengths of weak stellar photospheric features, 
and therefore direct systemic redshifts. Instead, galaxy 
systemic redshifts were esti r nated with the formulae pre- 
sented in lAdelberger et aP l)2003() . based on stronger in- 
terstellar features, and used to shift the spectra into the 
rest frame. 

3. RESULTS 

3.1. Direct Detection of Lyman- Continuum Emission 

Figures ^ a-nd [3 show rest-frame ultraviolet spectra 
for the sample over the rest-frame wavelength range 
800 — 1600 A. Blue and red side spectra have been 
joined at the dichroic wavelength, which ranges in rest- 
frame wavelength from 1175 A for the highest redshift 
object in the sample to 1345 A for the lowest redshift 
object. In each spectrum, the most prominent spec- 
tral features have been marked. The strongest feature 
is H I Lya, which appears in a variety of forms in our 
sample. The Lya profiles range from strong, complete 
absorption (Wo(Lya) ~ — 40A) to a mixture of emis- 
sion and absorption (IVo(Lya) ~ OA), to strong emis- 
sion (Wo (Lya) ^ 30A). The large range of Lya profiles 
and equivalent widths in z ^ 3 LBGs has been exam- 
ined previously (Shapl ev et al .112005) . and may provide 
important clues about the properties of the ISM in these 
galaxies. In an effort to showcase interstellar metal ab- 
sorption lines and the extremely faint Lyman-continuum 
region, while placing all of the objects on similar verti- 
cal scales, we truncated the full vertical scale of Lya for 
the six objects with the largest Lya emission equivalent 
widths (D3-apl, D17, MD23, C24, Cll, C32). In addi- 
tion to Lya, we mark low-ionization interstellar absorp- 
tion lines from Si II A1260, O I A1302 + Si II A1304, C II 
A1334, and Si II A1526, and high-ionization features from 
Si IV AA1393, 1402 and C IV AA1548, 1550. The Si IV 
and C IV features also contain contributions from stellar 
wind absorption and emission. Finally, the H I Lyman 
limit at 912 A is indicated. We note that two spectra are 
included for the object, D3 (D3-apl and D3-ap2). This 
object is characterized by two distinct components sep- 
arated by 200 km s~^ in velocity, and 1'.'9 on the sky 
15 proper kpc at z — 3.07). 
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Fic;. 1. Deep LRIS spectra of LBGs in the SSA22a field. The spectra have been shifted into the rest frame, and the flux units arc /t,h'. 
The strongest spectral features have been indicated in each panel, including Lya cinissioii/absorptioii, absorption from Ly/3 and higher-order 
Lyman series, and from low- and liigli-ionization interstellar lines. The Lyman limit is also marked. In this panel, the spectrum of D3-apl 
is the only one in which significant Lyman-continuum emission was detected. Typical exposure times in the Lyman-continuum region are 
8 hours using the UV-optimized LRIS-B setup, and the 1-cr uncertainty on the average flux density in the 880 — 910 A wavelength range 
(due to pixel-to-pixel fluctuations) is ~ 1 x lO"^'^ erg s~^cm~^Hz~^. 
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Fig. 2. — Deep LRIS spectra of LBGs in the SSA22a field (continued). Units, labels, and exposure times as in Figure [Tl In this panel, 
the spectrum of C49 is the only one in which significant Lyman- continuum emission was detected. 



Fig. 3. — Two-dimensional spectra of D3 and C19, objects with Lyman-continuum detections. Wavelength increases from left to right, 
and the spatial scale along each slit is indicated with a vertical bar of 5" in extent. The displayed portions of the spectra span in observed 
wavelength from 3350 — 3950 A, corresponding to ~ 830 — 970 A and ~ 810 — 950 A in the rest frames of D3 and C49, respectively. In 
each two-dimensional spectrum, the continuum clearly extends below the observed wavelength of the Lyman limit, which is indicated with 
an arrow. D3 actually consists of two components, as is apparent at the long-wavelength (right-hand) edge of this spectrum. However, the 
lower component (D3-ap2), exhibits no significant flux in the Lyman-continuum range. 
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Fig. 4. — Zoomed-in one-dimensional spectra of D3 and C49, objects with Lyman-continuum detections. The region extending down 
from Lya through the Lyman-continuum region is featured. Multiple higher-order Lyman-series are significantly detected in C49. The 
weakness of Ly/3 absorption in D3-apl is striking, however, as is the lack of a clear break below the Lyman limit. We note that Haro 11, 
the nearby starburst galaxy with a recent Lyman-continuum detection, is also characterized by the lack of a strong b reak at 912 A in th e 
rest-frame ( Bergvall et al...2005.V Motivated by the 2^3 average ionizing photon attenuation length of Az ~ 0.18 IMadau et aljlT99HI^ . 
we have generally evaluated the level of flux in the Lyman-continuum region over the narrow range of 880 — 910 A in this work. However, 
both D3-apl and C49 exhibit significant Lyman-continuum emission well below 880 A. 
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The measurement of escaping ionizing radiation con- 
stituted one of the principal objectives in obtaining such 
deep spectra of z ~ 3 LBGs. In order to measure the 
average flux density directly below the Lyman limit, 
fflnn, w e adopted the strategy presented in lSteidel et alJ 
l)2001t) . The rest- frame wavelength range used for evalu- 
ating /goo was 880 — 910 A, chosen to match theoretical 
estimates of the mean free path, Az 2± 0.18, for Lyman- 
continuum photons at z ~ 3 (Madau ct al. 1999). In gen- 
eral, the evaluation of /goo over a wavelength baseline sig- 
nificantly longer than the mean free path of ~ 30 — 40 A 
in the rest frame would be more indicative of the effects 
of opacity from intervening intergalactic systems than 
from the intrinsic ISM of the galaxies. 

Closer inspection of Figures ^ and |21 reveals signifi- 
cant Lyman-continuum flux in two spectra, those of D3- 
apl (z = 3.067), and C49 (z = 3.151). For D3-apl, 
we find /goo = 11.8 ± 1.1 x 10~^^ erg s-icm-2Hz"\ 
while for C49, we measure /goo — 6.9 ± 1.0 x 
10~^^ erg s~^cm~^Hz~^. The l-cr errors listed here cor- 
respond to statistical, pixel-to-pixel, fluctuations. In 
both of these objects, the detection of flux below the 
Lyman limit is signiflcant enough to stand out in two- 
dimensional, unextracted spectra, which are shown in 
Figure |3| A zoomed-in one-dimensional view of the 
Lyman-continuum regions is provided in Figure 0] In 
both the two-dimensional and zoomed-in one dimen- 
sional spectra, it is apparent that leaking flux exists at a 
roughly constant level well past 880 A, down to ~ 810 A 
in D3-apl, and ~ 800 A in C49. As mentioned above, 
the object D3 consists of two distinct continuum com- 
ponents. While the component D3-apl exhibits signif- 
icant flux at wavelengths shorter than 912 A, D3-ap2 
has no corresponding flux, indicating the complex na- 
ture of this source, and the limits on the physical size 
of the region where Lyman-continuum emission is escap- 
ing. In the 13 remaining spectra (12 additional objects, 
plus D3-ap2), there are no signiflcant detections of emis- 
sion below the Lyman limit. Among individual spec- 
tra in our sample, the average Icr uncertainty in /goo is 
fgoo, Stat = 1-15 X lO"'^-^ erg s~^cm~^Hz~^, which cor- 
responds to iriAB = 28.75, and a 3cr upper limit of 
ruAB — 27.56. These numbers, again, only represent 
pixel-to-pixel fluctuations. 

As shown in Figure upon closer inspection we flnd 
that the average flux density at 880 — 910 A for the ob- 
jects without signiflcant detections is distributed around 
a value that is slightly, but signiflcantly, negative. This 
artifact most likely reflects difficulty in achieving good 
near-UV fiatfields, free from systematic spatial fluctu- 
ations. Such systematic undulations are visible in the 
two-dimensional spectra of some of the objects with- 
out signiflcant Lyman-continuum detections. In order 
to correct for this systematic effect, we make the as- 
sumption that the objects without signiflcant detections 
should have Lyman-continuum flux densities distributed 
around zero, and add a constant value of A/goo = 
1.3 X 10"'^^ erg s~^cm~^IIz~^ to all spectra in the sam- 
ple. Accordingly, the Lyman-continuum flux densities 
listed above for C49 and D3-apl also include this sys- 
tematic correction. The offset applied represents a lower- 
limit, merely the amount required to make the non- 
detections distributed around zero ffux, rather than an 
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Fig. 5. — Sample distribution of flux densities in the Lyman- 
continuum wavelength range, 880-910 A. The two significant de- 
tections, for C49 and D3-apl, are clearly offset from the remainder 
of the sample. Objects without significant detections have flux 
densities distributed around an unphysical, negative value, which 
is significantly different from zero when considered as the average 
of 13 non-detections. To ensure that the objects without signif- 
icant detections are distributed around zero, we applied a sys- 
tematic additive correction to all of the spectra, at the level of 
A/900 = 1-3 X 10~^^ erg s~lcm~^Hz~^. As the objects without 
significant detections may also be distributed around a value that is 
slightly positive, the correction we applied represents a lower limit. 
The la width of the distribution for objects without significant 
detections indicates the total level of uncertainty on an individual 
measurement of /900, including contributions from both pixel to 
pixel and systematic errors. 



unphysical, negative value. We cannot rule out the pos- 
sibility that the individual non-detections are distributed 
around a small, slightly positive value, and therefore re- 
quire a larger correction, but for the remainder of the 
discussion we will adopt the conservative lower-limit to 
interpret the data. To place the described correction in 
context, we note that it represents on average 2% of the 
continuum level at 1500 A for objects in our sample. 

The width of the distribution of Lyman-continuum flux 
densities shown in Figure indicates the total level of 
uncertainty in the individual /goo values, including both 
pixel-to-pixel uncertainties and systematic sources of er- 
ror arising from sky subtraction and flat-flelding. This la 
width is 1.9 x 10^'^^ erg s^^cm^^Hz^"'^. Given the typical 
value of the /goo pixel-to-pixel uncertainty, trgoo, stat, we 
flnd that the contribution to the error from systematic 
sources is crgoo.sys 1-5 x 10"'^^ erg s~^cm~^Hz~^. This 
value pertains to /goo values measured from individual 
spectra and must be divided by the square-root of the 
number of spectra for the associated composite spectra. 
Throughout the remainder of the current work, we sep- 
arately indicate statistical pixel-to-pixel, and systematic 
flat-flelding- and sky-subtraction-related errors. 

As well as considering the Lyman-continuum flux from 
individual objects in our deep sample, we also measured 
the properties of a composite spectrum constructed from 
the whole sample, shown in Figure This spectrum 
represents a simple average of all of the individual spec- 
tra following their systematic correction as described 
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Fig. 6. — Composite spectrum of deep LBG spectra. This plot represents the average of the sample of spectra displayed in Figures 111 
and 13 and shows flux in the Lyman-continuum region that is formally significant at greater than a 3cr level. The sample average ratio of 
1500 A to Lyman-continuum flux density, uncorrected for IGM opacity, is (/1500//900) = 58 it ISstat ± 17 sys- 



above. Given the manner in which the spectra were 
corrected, the mean /goo represents the contributions 
from C49 and D3-apl alone, averaged over the whole 
sample. In the average composite spectrum, we find 
/900 = 1.08±0.28,tat±0.37,j,,x 10-31 erg s-icm-^Hz-^. 
While it does not reflect the total uncertainty in the 
Lyman-continuum flux-density level, the pixel to pixel 
noise represents a Icr depth ofruAB = 30.3, almost a fac- 
tor of 3 deeper than the corresponding itiab = 29.1 la 
depth of the composite LBG spectrum used to m easure 
Lyman-continuum emission in lSteidel et alJ l|20fllD . 

3.2. Properties of Galaxies with Lyman- Continuum 
Detections 

We now consider the observed properties of D3 and 
C49 in order to discern if they stand out from ob- 
jects without significant Lyman-continuum flux detec- 
tions. There are a number of additional observations of 
D3, including ground-based and HST/NICMOS near-IR 
imaging, ground-based near-IR spectroscopy, and deep 
narrow-band imaging tuned to the wavelength of Lya 
at the redshift of the SSA22a galaxy overdensity. There 
were no spectroscopic observations of C49 prior to the 
deep ones presented here, and there are accordingly fewer 
multi- wavelength data points for this object. FigureElin- 
cludes plots of several different rest-frame UV properties 
for objects in the deep sample, including Lya equiva- 
lent width; average low-ionization interstellar absorption 
equivalent width (based on the strengths of Si II A 1260, 
O I A1302 + Si 11 A1304, C II A1334, and Si 11 A1526); 
E{B—V)] and rest-frame UV absolute magnitude, uncor- 
rected for dust extinction. D3 and C49 are highlighted 
in each figure, as well as the average values for the three 



different LBG samples listed in Table 2. 

3.2.1. Rest-frame UV Properties 

D3 has the brightest rest-frame UV luminosity in the 
current sample, with TZab — 23.37. The flux in the TZ- 
band represents the sum of the flux from both compo- 
nents of D3 (i.e. D3-apl and D3-ap2), which were treated 
as a single object in the original photometric catalog. As 
estimated from the ratio of continuum levels in the spec- 
tra of D3-apl and D3-ap2, D3-apl is approximately two 
times as bright as D3-ap2, and this ratio is roughly con- 
stant over the wavelength range 4000-6300 A, where the 
flux-calibration is accurate. The ratio between the D3- 
apl and D3-ap2 spectra should closely resemble the true 
component ratio, since the slit position angle of —50° 
east of north is within 10° of the observed angle between 
the two components in the optical imaging data. C49 
is brighter than the median of the deep sample in rest- 
frame UV luminosity, while not as extreme as D3. In 
a broader comparison of rest-frame UV luminosity, D3 
is among the brightest 2% of the 811 obj ects in the to- 
tal LBG spectroscopic sample presented in lShaolev et alJ 
||^03J); C49 would have been among the brightest 8%. 

We have interpreted the rest-frame UV colors of galax- 
ies in the current deep sample assuming that the unred- 
dened continuum is well-represented by a 300 Myr con- 
tinuous star-formation model, where 300 Myr repre- 
sents the med i an ag e in the LBG sample analyzed by 
iShaolev et al.l l)2001li . Differences in G — 7?. color (cor- 
rected for IGM absorption and Lya emission/absorption 
intrinsic to the galaxies) are interpreted as differ- 
ent amounts o f dus t extinction, parameterized by a 
iCalzetti et al.l l)2000() attenuation curve. Accordingly, 
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Fig. 7. — Rest-frame UV properties of the deep spectroscopic sample. In each panel, solid black points indicate individual objects 
in the deep spectroscopic sample; points for C49 and D3 are enclosed in open squares; the large open circle, star, and triangle indicate 
the average values for the deep spectroscopic sample, the total LBG spectroscopic sam ple from lS hanlev et al.l 120031) . and the original 
Lyman-continuum sample of Steidcl ct al. (2001), respectively. In general, while the Stcid el et alj 120011) sample is very similar to the total 
LBG sample in terms of average Lya and low-ionization interstellar equivalent widths, the current sample is much closer to the total LBG 
sample in terms of average UV extinction properties. (Upper left) The correlation between total Ly« equivalent width and low-ionization 
interstellar absorption equivalent width in composite spectra (Shapley et al.,|.20p3!) is confirmed on an object-by-object basis for the deep 
sample. Positive values of equivalent width indicate emission. (Upper right) The correlation between total Lyo equivalent width and 
UV extinction. D3 has the largest inferred value of E{B — V) in the deep sample. (Lower left) The correlation between E{B — V) and 
low-ionization interstellar absorption equivalent width. In this plot, D3 stands out for having an anomalously high dust extinction, given 
the equivalent width of its low-ionization interstellar absorption lines. (Lower right) A plot of total Lyo equivalent width as a function of 
UV absolute magnitude , uncorrected for dust extinction. D3 is the brightest object in the sample. Also, it is clear that both the current 
and lSteidel et aLI 120011) samples are significantly more luminous on average than the total LBG spectroscopic sample. 



the observed G—TZ colors can be converted into E{B — V) 
values, with the largest associated uncertainties resulting 
from the uncertainties in assumed underlying stellar pop- 
ulation template, and the observational errors on G — TZ 
color (typically at the 5 — 10% level). We note that es- 
timates of UV-continuum slope and E{B — V) inferred 
directly from the shapes of the deep spectra redwards of 
Lya are not very robust, because the changes in relative 
flux level corresponding to different degrees of UV ex- 
tinction are not significant over the range covered by the 
spectra. 

As shown in Figure[7| D3 has E{B -V) ^ 0.23, which 
implies the highest level of extinction in the deep sam- 
ple, more than 2 magnitudes of attenuation at rest-frame 
UV wavelengths. D3 is significantly redder than the av- 
erage in the deep sa mple, the total LBG sample from 
Sliagle^_et_^lJ ('2003') , and the sample of LBGs used by 
Steidel et ^(2001) to construct a composite spectrum 
that showed significant Lyman-continuum emission (see 
Table 2). In fact, the level of dust extinction in D3 may 
be even higher; simu ltaneously fitting th e G ~ TZ and 
n-K, colors of D3 l|ShaDlev et al.ll200H) yields an age 



< 10 Myr, and E{B — V^) > 0.3 (i.e. the assumption of 
300 Myr for the age will lead to an underestimate of the 
extinction if the true age is younger). On the other hand, 
Redd v et all lj200(y find that LBGs with the youngest 
inferred ages do not appear to follow the iMeurer et alJ 
(199!|) relation between UV-reddening and dust obscu- 
ration, in the sense that they have redder UV colors for 
a given amount dust extinction. The young apparent 
age of D3 might therefore cause us to overestimate the 
true level of dust obscuration. As we return to in sec- 
tion^ the heavy dust extinction inferred from the UV 
colors of this object is somewhat difficult to reconcile 
with the sign i ficant Lyman-continuum emission. In fact, 
'Steidel et al.l 12001^ propose that it is galaxies with less 
reddened continua that are more likely to show escap- 
ing Lyman-continuum emission. We have also checked 
whether the composite nature of D3 could have led to a 
bias in the measured G — TZ color, but find that the mor- 
phology between the G and TZ bands remains roughly 
constant (which is consistent with the roughly constant 
ratio between the continua of D3-apl and D3-ap2 spectra 
spanning from 4000 A to 6300 A). Therefore, the 7e-band 
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detection isophote applied to the G band should yield an 
accurate representation of the G — TZ color. In contrast to 
D3, the other object with significant Lyman-continuum 
emission, C49, has an inferred E{B — V) = 0.06. This 
low value of E{B — V) places C49 among the third of the 
deep sample with the least extinction. It appears signifi- 
cantly less reddened than the average of the deep sa mple 
and the total LBG sample of iShaplev et al.l l)2003D . and 
even slightly bluer than the aver age of the blue sample 
presented in iSteidel et all l)200H) (see Figure [71). Thus, 
no consistent picture emerges from D3 and C49 about 
the dust extinction properties of galaxies with significant 
escaping Lyman-continuum emission. 

Now we compare the rest-frame UV spectroscopic fea- 
tures of D3 and C49 with those of the samples listed 
in Table 2. Lya is observed in emission for both com- 
ponents of D3, with a rest-frame equivalent width of 
~ 5 A in both apertures. C49 also shows weak Lya 
emission, but superposed on an absorption trough of 
roughly equal strength, such that the total equivalent 
width is ^ OA. The Lya emission in C49 is actually 
double-peaked, though this fact is not readily apparent in 
Figure 01 which has been smoothed by a 3-pixel boxcar. 
As shown in Figure [7| D3 and C49 exhibit close to the 
average Lya equivalent width in the deep sample, which 
is itself characterized by weake r Lya emission on ave rage 
than the larger LBG sample of lShaolev et alJ l|2003f l and 
the sample used to detect Lyman-continuum emission in 
ISteidel et a,1. (2001). 

Another relevant property is the strength of the 
low-ionization interstellar metal absorption lines. The 
strongest low-ionization interstellar absorption lines in 
the spectrum of D3 (aperture 1) are consistent with zero 
intensity at line center; in this aspect, the D3 spectrum 
may imply a unity covering fraction of gas giving rise to 
low-ionization metal absorption, at least at some veloci- 
ties. As shown in Figured the low-ionization absorption 
equivalent widths in D3 are close to the mean of the deep 
sample, which is stro nger than the ave r age in both the 
IShaplev et all l|2003 D and I Steidel et all l|2001f ) samples; 
D3 also has narrower low-ionization velocity widths than 
the deep sample average, as estimated from the Si II 
A1260 and C II A1334 absorption lines. Also shown in 
Figured C49 has weaker than average low-ionization ab- 
sorption lines that do not reach zero intensity at line cen- 
ter, even after deconvolution of the spectral point spread 
function. In the absence of unresolved saturated com- 
ponents for which we cannot correct, it appears that 
C49 is characterized by neutral gas with non-unity cov- 
ering fraction - i.e. a patchy spatial distribution. The 
strength of the low-ionization interstellar absorption lines 
in the spectrum of C49 is very simila r to those m easured 
from the composite spectrum of Steidel et alJ 12001.) , 
in which significant Lyman-continuum emission was de- 
tected. Weak low-ionization interstellar absorption does 
not however appear to be a sufficient condition for escap- 
ing Lyman-continuum emission. Three galaxies in the 
deep sample (C32, D17, and MD23) have weaker low- 
ionization absorption lines than C49, and C32 and D17 
also appear to have lower covering fractions of neutral 
gas. However, none of these three galaxies exhibit any 
significant emission in the Lyman-continuum region. 

One striking property of D3 is the apparent weakness 
of Lyman series absorption. The noise in the spectrum 



at these wavelengths prevents a precise determination of 
H I column density, yet a comparison with the average 
spectrum of the deep sample reveals significantly weaker 
than average absorption in D3 from Ly/3, Ly7, Lye. In 
particular, the component of D3 with detected Lyman- 
continuum emission is quite distinct from the rest of the 
deep sample in terms of the weakness of Ly/3 absorption. 
On the other hand, C49 agrees very well with the deep 
sample average for Ly/3 and the higher-order Lyman se- 
ries absorption features. 

The spectra of D3 and C49 feature both Lya emission 
and interstellar absorption lines, and the higher redshift 
of Lya emission relative to that of the absorption lines 
indicates kinematic evidence for large-scale outflow mo- 
tions. This kinematic signature appears to be an almost 
universal char acteristic of LBG re st-frame UV spectra 
(jPettini et all [2001: Sh aplev et all 12003). The velocity 
offset between Lya emission and interstellar absorption 
lines is ^ 500 km in D3 (aperture 1) and some- 
where between 780 and 1130 km s~-^ in C49, where the 
range for C49 indicates the different offsets calculated 
using each of the two peaks of Lya emission. The veloc- 
ity offset in D3 is smaller than the deep sample average 
for objects with both emission a nd absorption, as we ll 
as the average offset measured bv 'Shapl ev et all ( )2003D . 
C49, on the other hand, shows a larger than average 
offset between emission and absorption velocities. Addi- 
tional rest-fr ame optical nebular emission line measure- 
ments of D3 ijPettini et al.ll2f)f)lf) indicate that the low- 
ionization interstellar absorption lines are blueshifted by 
~ 280 km s^^ relative to the H II regions, while Lya 
is redshifted by only ~ 215 km s~^. Therefore, in D3, 
it is the redshift of Lya relative to H II regions which 
is smaller than the typical value found in LBGs, not 
the blucshift of the low-ionization interstellar absorp- 
tion lines, and the associated outflow speed (Pcttini et aO 
2001: Steidel et al.ll200l . It wiU be valuable to obtain 
H II region emission line measurements for C49 as well, 
in order to locate the systemic redshift for this system, 
and determine the outflow speed from the blueshift of 
interstellar absorption lines. 

Although the higher-ionization interstellar absorption 
lines are noisy, we can roughly estimate that the Si IV 
AA1393, 1402 absorption is comparable in D3 to that in 
the average spectrum of the deep sample. The Si IV 
profiles are too noisy to determine a useful estimate of 
the velocity widths. It is not possible to measure the 
strength or width of the C IV line in D3, as the feature 
is contaminated by a strong sky residual. In C49, both 
Si IV and C IV interstellar absorption are weaker than 
the average of the deep sample. Additionally, the Si IV 
lines are significantly narrower than the average Si IV 
velocity width in the deep sample. 

3.2.2. Additional Properties of D3 

D3 is one of the brightest objects in the total spec- 
troscopic LBG sample, and, as such, has served as an 
obvious target for multi-wavelength follow-up. The com- 
bination of Keck/NIRC i^g-band photometry with exist- 
ing optical observations indicates a young age for D3, of 
less than 10 Myr, and a s tellar mass of ~ 6 x 1O^M0, 
assuming a Chabricr IMF f Shaplev et"an i20011. The in- 
trinsic ratio of Lyman-continuum to non-ionizing UV lu- 
minosity is a decreasing function of age, so that the high- 
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Fig. 8.— n (left) and i?5T/NICMOS F160W (right) images of D3. The TL image was obtained at the Palomar 200 inch Hale Telescope, 
under ~ 1" seeing conditions, while the F160W image was obtained with the NIC2 camera of NICMOS, which provides diffraction-limited 
(0'.'14 resolution) imaging at 1.6/im. The locations of spectral apertures 1 and 2 are marked in each image, as "apl" and "ap 2," respectively. 
The vertical bar in the lower left-hand corner of each figure indicates a scale of 5" . D3-apl is the component for which flux was detected in 
the Lyman continuum region. Both the TZ. and H5T/NICMOS F160W images indicate the brighter and more compact nature of D3-apl, 
relative to D3-ap2. D3 is the only object in the deep sample with confirmed double morphology. 



est intrinsic ratio would be found in galaxies with ages 
as young as that inferred for D3 - though, of course, 
the true age for D3 cannot be less than 10 Myr, i.e., 
shorter than its dynamical timescale. Observations of 
C49 in the near-IR constitute a natural follow-up pro- 
gram and will allow similar age and stellar mass deter- 
minations. Keck/NIRSPEC Near-IR spectroscopy cover- 
ing the [OIII] AA5007, 4959 H/3, and [Oil] A3727 emission 
lines enab les an estima te of the -R23 chemical abundance 
indicator llPagel et alJll979ll . Based on the value of i?23 , 
IPettini et abl lj2nnir> find O/H=0.24 - 0.78 (0/H)q for 
D3, where these values refer to the rev ised solar oxy- 

fQ abundance of 12-|-log(0/II)Q — 8.66 (jAsplund et al.l 
04). The NIRSPEC spectrum is roughly centered on 
the D3 component showing significant Lyman-continuum 
emission, and the slit was not positioned at an angle to 
cover both components. Further near-IR information in- 
cludes three orbits of /f^T/NICMOS F160W diffraction- 
limited imaging. As shown in Figurc|Sl the rest-frame op- 
tical image confirms with significantly sharper resolution 
what is suggested in the TZ band image of D3. Specif- 
ically, the component with escaping Lyman-continuum 
emission has a peak surface brightness 2.5 — 3.0 times 
higher than that of the other component. The brighter 
component is also more compact and radially-symmetric, 
with a diameter roughly half the extent of the longest 
linear dimension of the extended, irregular and diffuse 
emission of the second component. Longer wavelength 
observations of D3 thus indicate the complex morphol- 
ogy of a young system with a warm ionized gas phase 
that has already been enriched to a significant fraction 
of solar metallicity. 

As the redshift of D3 places it within the well-studied 
overdensi ty of LBGs at z = 3.09 ± 0.03 (Stcidcl et al. 
[1,998, 2000), narrow-band imaging can be used to reveal 
the spatial distribution of Lya e mission in this object . 
The narrow-band observations of iMatsuda et al.l l)2004fl 
indicate that Lya emission is extended over 17 arcmin^ 
with a luminosity of 9.5 x 10'*^ erg s~^, and roughly cen- 
tered on the 7^-band position of the component with de- 
tected Lyman-continuum emission. The large spatial ex- 
tent of its Lya emission places D3 in the so-called Lya 



"blob" sample of IMatsuda et all lj2004(l . The Lya lu- 
minosity reflects the number ionizing photons absorbed, 
and re-emitted as Lya photons by recombining gas. For 
more than one third of the "blob" sample, the num- 
ber of intrinsic Lyman-continuum photons inferred from 
the UV continuum is exceeded by the number inferred 
from the Lya luminosity. Potential explanations for 
the excess of Lya emission include photoionization by a 
young, metal-poor stellar population with a flatter than 
Salpeter IMF at the high-mass end; by an AGN; by hid- 
den ionizing UV sources; or by the dif fuse intergalac- 
tic ba ckground radiation. Figure 10 of IMatsuda et alJ 
lj2004(l . showing the star-formation rate inferred from 
Lya compared with that inferred from the non-ionizing 
UV-continuum, indicates that no such explanation is re- 
quired for D3. The number of ionizing photons repro- 
cessed into Lya emission represents only < 25% of the 
number inferred from the non-ionizing UV continuum. 
Understanding the extended Lya emission in D3 is yet 
another goal in untangling the complex nature of this ob- 
ject. Since the redshift of C49 does not place it within the 
SSA22a galaxy overdensity, no information is presently 
available about the spatial distribution of its Lya emis- 
sion. 

In summary, no consistent picture emerges from D3 
and C49 about the properties of galaxies with escap- 
ing Lyman-continuum emission. While C49 is one the 
bluest galaxies in the deep sample, D3 has the reddest 
observed UV continuum. While the low-ionization inter- 
stellar metal absorption lines from C49 indicate a non- 
unity covering fraction of cool gas, favorable for allowing 
the escape of Lyman-continuum photons, there are galax- 
ies in the deep sample with weaker low-ionization inter- 
stellar absorption lines and also apparently lower cover- 
ing fractions of cool gas, which at the same time do not 
exhibit significant escaping ionizing radiation. Further- 
more, D3, which does exhibit escaping Lyman-continuum 
radiation, has low-ionization interstellar absorption lines 
that may be black at line center, indicating a unity cov- 
ering fraction of cool gas. Two striking properties of 
D3 are its complex double morphology, consisting of a 
compact component exhibiting Lyman-continuum, and a 
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fainter and more diffuse component with no correspond- 
ing detection; and its weak Lyman series absorption - 
significantly weaker than the typical level of absorption 
in the deep sample. At the same time, C49 has a sim- 
ple single-component morphology, and Lyman series ab- 
sorption that agrees quite well with that in the sample 
average spectrum. Clearly, a sample of two galaxies is 
too small to infer the typical properties of galaxies with 
escaping Lyman-continuum emission. 

4. THE ESCAPE FRACTION OF IONIZING RADIATION 

The detection of escaping Lyman-continuum emission 
from individual star-forming galaxies at z ~ 3 represents 
definite forward progress from previous statistical detec- 
tions (Steidel et al. 2001) a nd indiv idual non -detections 
itGiallongo et al... .2002: Inoue et alJ 12005.: .Dawson et all 
120021) . The range of detections and non-detections of 
ionizing flux in our sample can also be combined with 
measured 1500 A fluxes to characterize the emergent far- 
UV spectral shapes of star- forming galaxies at z ~ 3. 
Ultimately, our goal is to quantify the global contribu- 
tion of LBGs to the ionizing background at z ~ 3. The 
most realistic route for addressing this problem is to con- 
vert the well-characterized LEG 1500 A luminosity func- 
tion ( Steidel et al. 1999t. lAdelberger fc Steidel' 20001 into 
a Lyman-continuum luminosity function using the aver- 
age emergent ratio of 1500 A to Lyman-continuum flux 
densities, corrected for IGM absorption, and determined 
from a small but representative s ample of LBGs. 

From a theoretical perspective (Loeb fc Bar kanall2001t 
iFurlanetto et all l2004t iCen fc Ostrikcr 1993 }), it is of 
interest to determine the absolute escape fraction of 
Lyman-continuum photons, fesc- This quantity, defined 
as the ratio of 900 A photons produced that escape from 
the galaxy without being absorbed, relative to the total 
number produced, is difficult to estimate from even the 
most sensitive observations, because it relies on knowl- 
edge of the intrinsic number of Lyman-continuum pho- 
tons pro duced. A re l ated q uantity is fesc,rei, first de- 
fined by Steid el et alJ (j2001) as the fraction of escaping 
Lyman-continuum photons normalized by the fraction of 
escaping 1500 A photons. In practice, fesc.rei is a better 
defined observational quantity, since it can be combined 
with the well-determined UV continuum luminosity dis- 
tribution to estimate the global contribution of LBGs to 
the ionizing background. It is also possible to determine 
fesc from fesc,rei if the degree of dust extinction in the 
UV is known. The relative escape fraction, fesc,rei, can 
be expressed in terms of the observed 1500 A to 900 A 
flux-density ratio, /1500//900, as: 

J. (il50o/-^90o)mt / s 

Jesc.rel — — 77 77 ^ — ^^P(TlG M, 900 ) Uj 

U1500//900j 

where (iisoo/^ooo)*™* is the intrinsic ratio of non- 
ionizing to ionizing specific intensities, and tigai,9oo rep- 
resents the line-of-sight opacity of the IGM to Lyman- 
continuu m photons, which ca n be estimated either em- 
pirically (|^ei^^^I]_|22^) or t hrough simu lations 
ijMadaul 119951 iBershadv etal]IT999t llnoue et alJl2005l) . 
The real difficulty in estimating fesc,rei and fesc stems 
from the fact that the intrinsic spectral shape for mas- 
sive stars above and below the Lyman Break is not 
well-constrained by observations. (ii5oo/-^90o)int must 



therefore be estimated from stellar population synthe- 
sis models, and varies in the models as a function of 
stellar population age, metallicity, star-formatio n his- 
tory, and IMF (Bruzual & Chariot 2003; Lcithcrcr et alJ 
11999! ). Even wh en fixing the IMF and star-formation his- 
tory, Qnou^OL (,2005.) find that (Li5oo/i90o)mt varies 
from 1.5 — 5.5, mainly as a function of time since the 
onset of continuous star-formation (where the units of 
- Lisnn and Lgnn are erg s~^H z~^). Both iSteidel et afl 
()2001|) and llnoue et al.l l)2005j) adopt an intrinsic ratio 
of (ii5oo/^90o)mt = 3.0. For comparison with previous 
work, we adopt the same ratio here, but note the signif- 
icant associated uncertainties. 

In this section, we present estimates of /1500//900 
and fesc,rei for the two individual detections, as well 
as an estimate of the sample average /1500//900 and 
fesc,rei ■ Modeling and correcting for the effects of IGM 
Lyman-continuum opacity enables the conversion of a 
raw /1500//900 value to the emergent /i5oo//90o,corr that 
would be observed in the immediate vicinity of a galaxy, 
which is one of the quantities required for estimating 
the global contribution of LBGs to the ionizing back- 
ground at z ~ 3. In order to model IGM opacity, we 
use Monte Carlo simulations of intergalactic absorption 
to generate large samples of random sightlines to the 
redshifts of C49 and D3, and to the average redshift of 
the deep sample. Absorbers are drawn at random from 
the column-density and redshift distributions constitut- 
ing the MC - NH model of ( Bershadv ct al. 1999), and 
then applied to the unabsorbcd continuum of each simu- 
lated sightline. The resulting sample of model sightlines 
is analyzed to obtain the mean and standard deviation 
of the optical depth in the rest-frame region 880 — 910 A 
at each emission redshift. As an example, we show the 
simulated distribution of optical depth for the sample av- 
erage redshift, (z) = 3.06, in FigurelHl In order t o corr ect 
the composite LBG spectrum in Steidel ct al. l)2001fl . a 
more empirical approach to estimate the IGM opacity 
at z = 3.47 was adopted. A sample of 15 bright QSO 
spectra at (z) = 3.47 were combined to create a compos- 
ite spectrum with no intrinsic Lyman-continuum opacity, 
yielding the average intervening absorption in the Lyman 
continuum over 15 random sightlines. The absorption es- 
timated in this empirical fashion was a factor of ^ 3.8, 
which agrees well with the estimate of a factor of ~ 3.7 
from our simulated z = 3.47 sightlines. Finally, in this 
section, we also attempt to quantify the sample varia- 
tion in /i5oo//900jcorr, taking into account sightline to 
sightline variations in IGM opacity. 

4.1. The Emergent UV Spectrum: C49, D3, and the 
Sample Average 

The observed UV-to-Lyman-continuum flux-density 
ratio in C49 is /1500//900 — 12.7 ± 1.8. Simulations of 
IGM neutral hydrogen opacity indicate that, on average, 
the flux level in the rest-frame range 880 — 910 A is at- 
tenuated by a factor of 2.8 for sources at z = 3.15 (the 
redshift of C49). If the hue of sight to C49 is typical 
for its redshift, then the emergent far-UV flux-density 
ratio from the galaxy is /1500//900, corr = 4.5, very sim- 
ilar to that i n ferred from the composite LBG spectrum 
ISteidel et alJ l)200H) . Furthermore, the rest-frame UV 
colors of C49 imply that differential dust extinction be- 
tween the UV and Lyman-continuum is not significant. 
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Fig. 9. — Results of Monte Carlo simulations of IGM opacity 
along the line of sight to z = 3.06, the average redshift of the deep 
sample. A large number N = 5000 of random sightlincs to 2 = 3.06 
were generated. For each sightline, the ratio of the attenuated 
to unattenuated flux densities in the rest-frame 880 — 910 A was 
computed, and is expressed here as exp(— rgoo). The attenuation 
factor discussed in the text consists of the inverse of exp(— rgoo). 
At z = 3.06, we find (exp(— rgoo)) = 0.38, corresponding to an 
attenuation factor of 2.6. The 95.4% confidence interval for the 
880—910 A attenuation factor at z = 3.06 spans fro m 35 to 1.6. The 
dotted line indicates the analytic approximation of IMadaul J1995I) 
for the mean transmission in the 880 — 910 A region at z = 3.06. 
The tail in the distribution extending down to zero transmission is 
dominated by sightlines in which at least one Lyman limit system 
fell within the redshift range corresponding to rest-frame 880 — 
910 A. 



Assuming th e same intrinsic ratio of (ii5oo/-^90o)int as 
ISteidel et alJ ll200l . we find f^sc,rei = 65% for C49. 
If C49 is, however, an unusually transparent sightline, 
the emergent UV-to-Lyman-continuum ratio is larger. 
Our Monte Carlo simulations of IGM opacity indicate 
that, at a 95.4% (i.e. 2a) confidence level, the at- 
tenuation in the 880 — 910 A rest frame falls between 
a factor of 38 and 1.7. At the transparent edge of 
the 2a confidence interval, the IGM attenuation implies 
/1500//900, corr= 7.5, and fesc,rei ~ 40%. Given that 
continuum flux is detected in C49 down to ~ 8OOA, it 
is likely for the transparency to be higher than aver- 
age in the ~ 880 — 910 A range. It is important to 
note here that one cannot infer the transparency level 
in the range ^ 880 — 910 A, based on Da: the continuum 
decrement at 1050 — 1170 A. Da indicates absorption 
from intervening Lya lines, and is dominated by low- 
column-density systems (i.e. N{HI) < 10^^ cm~^). On 
the other hand, the Lyman-continuum attenuation fac- 
tor is modulated by the number of intervening systems 
with N{HI) > 10^'' cm~^ at redshifts such that the as- 
sociated Lyman-limit absorption falls in the window of 
interest. Accordingly, our Monte Carlo IGM simulations 
indicate the absence of a correlation between Da and the 
attenuation factor at ~ 880 — 910 A. 

In D3, the UV-to-Lyman-continuum flux-density ratio 
is estimated using the 1500 A flux from both components, 
compared with the Lyman-continuum detection from 
aperture 1, alone. This comparison yields /1500//900 = 



7.5 ± 1.0. Correcting for the average 880 - 910 A at- 
tenuation of a factor of 2.6 at z = 3.07 (the redshift 
of D3), we find /1500//900, corr = 2.9, which is roughly 
the same as the intr insic ratio of L 150 /^909 assumed by 
ISteidel et ail l|200lD and llnoue et~aL[ f2005). and there- 
fore implies fesc.rei > 100%! Again, it is worth noting 
that the inferred significant dust extinction in D3 is dif- 
ficult to reconcile with such an extreme value of fesc.rei, 
as a reasonable extrapolation of the[Calzetti et al. ( 200(4) 
attenuation law would imply significantly more attenua- 
tion at 900 A than at 1500 A. 

There are multiple ways to explain the extreme and 
apparently contradictory properties of D3. First, in ad- 
dition to the possibility that interpreting the UV color 
with the Calzetti law causes us to overestimate the over- 
all dust extinction in D3 (see section r^.2.1|l . the E{B~V) 
inferred for this object may only suffice as a global quan- 
tity, averaged over the entire surface area of the galaxy. 
If, in fact, the escaping Lyman-continuum emission is 
patchy, and preferentially escaping from regions with less 
dust extinction, it may not be appropriate to use the av- 
erage E{B — V) to interpret the emergent ratio of UV to 
Lyman-continuum radiation. Future HST/ ACS imaging 
in two filters may constrain the detailed dust extinction 
distribution in this source. Second, the average correc- 
tion for IGM opacity along the line of sight may not 
apply to an individual sightline. For a given emission 
redshift, there is significant variation in the rest-frame 
880 — 910 A opacity from sightline to sightline. In our 
Monte Carlo simulations of IGM opacity, we found that, 
at 95.4% (2(7) confidence, the attenuation in the Lyman 
continuum at 2: = 3.07 ranged between 35 and 1.6. Us- 
ing 1.6 as a lower limit on the IGM attenuation, we 
derive fi5oo/ f 900, corr — 4.7. Compared to an assumed 
(Li50o/i^90o)mt ratio of 3, this yields fesc,rei ^ 65%. 
As in the case of C49, adopting a lower than average 
IGM opacity in the Lyman continuum may indeed be 
suggested by the persistence of the D3 continuum down 
to rest wavelengths of ~ 810 A, well below one Lyman- 
continuum mean free path. Third, we consider the possi- 
bility that the intrinsic UV continuum of D3 is not dom- 
inated by stars, but by the emission from an AGN, with 
a smaller intrinsic ratio of (^1500/^900 )int- This scenario 
is not supported by rest-frame UV spectroscopic obser- 
vations. We detect no strong C IV or N V emission in 
the rest-frame UV spectrum, which would be characteris- 
tic. Furthermore, HST/NICMOS high-resolution imag- 
ing reveals no unresolved nuclear emission in the rest- 
frame optical. Future X-ray and mid-IR observations 
with Chandra and Spitzer, respectively, may constrain 
the contribution of an AGN to the emission in D3. 

Finally, there is the possibility that the apparent de- 
tection of fiux below the Lyman limit in D3 results from 
scattered light from the outer isophotes one or both of 
the bright galaxies within several arcseconds of the posi- 
tion of D3-apl. These bright neighbors stand out in the 
TZ and F160W images in Figure|Hl with centroids 5'.'4 to 
the northeast of D3-apl, and 5'.'9 to the southwest. Per- 
haps most relevant, the brighter southwest object has 
Un = 21.6, while the northeast object has C/„ = 23.5. 
D3 is unique among the objects in the deep sample, in 
terms of being flanked by such bright objects. While 
the mask position angle was chosen to be roughly paral- 
lactic for exposures at hour angles east of overhead, the 



15 



direction of differential refraction for exposures obtained 
when the SSA22a field was setting is such that light from 
the brighter southwest object was bent in the direction 
of the D3 slit. However, even in the worst case scenario, 
the deflection at 3600 A for the setting exposures with 
the maximum airmass is still only 0'.'64, a small fraction 
of the distance between the bright southwest object and 
the edge of the D3 slit. Nonetheless, the unusual config- 
uration of D3 with its bright neighbors should be noted. 

We also consider the mean observed UV-to-Lyman- 
continuum flux-density ratio for the entire sample of 14 
objects with deep spectra, by constructing a composite 
spectrum for which each individual galaxy spectrum was 
scaled to a common mode at 1400 — 1500 A before averag- 
ing. This mean ratio is (/1500//900) = 58±18stat±17sys, 
where the two sources of error reflect the pixel to pixel 
noise in the composite spectrum and the systematic un- 
certainties in the individual values of /goo due to flatfield- 
ing and sky-subtraction, respectively. The mean red- 
shift of the sample is {z) = 3.06, very close to that 
of D3. Therefore, the same correction factor for IGM 
opacity (2.6) can be applied to this ratio to obtain the 
average emergent UV-to-Lyman-continuum ratio. With 
this correction, we obtain (/1500//900 icorr } — 22. As- 
suming an intrinsic UV-to-Lyman continuum ratio of 
(^i50o/-^90o)int — 3.0, as in ISteidcl ct al. (2001), we find 
{fesc,rei) = 14% for the deep sample - a f actor of ^ 4.5 
lower than the value of fesc,rei reported in lSteidel et al.l 
(jlOOl). We will return to this result in sectional While 
significantly smaller than the uncertainty in opacity for 
an individual sightline, the average opacity for a sam- 
ple of 14 sightlines has an associated error, because of 
the finite number of sightlines being averaged together. 
To estimate this uncertainty, we drew a large number of 
random samples of 14 sightlines from our simulations of 
single IGM sightline opacities. For each random sample, 
we computed the average IGM opacity, and then com- 
puted the mean and width of the distribution of IGM 
opacities averaged over 14 sightlines. With 95.4% {2a) 
confidence, we found that the average IGM opacity for 
the sample of 14 sightlines requires a correction factor 
between 2.2 and 3.3, implying that (/1500//900 icorr ) falls 
between 17 and 27. The corresponding range in fesc,rei 
is 11%-18%, assuming (ii5oo/i9oo)mt as before. 

There are two important caveats associated with 
our attempts to model and correct for IGM opacity 
and infer the shape of the emergent far-UV spectrum, 
fi50o/ f 900, corr- First, the IGM simulations we performed 
to characterize the distribution of attenuation factors in 
the Lyman-continuum region are based on the assump- 
tion that each random sightline is a probing the aver- 
age properties of the IGM, and do not include the ef- 
fects of large-scale structure. In fact, nine of 14 of the 
galaxies in our sample have redshifts that place them 
within the significant ovcrdcnsity in the SSA22a field 
discovered by Stcidcl ct al. (1998), which is character- 
ized by Sgai - 6 at (z) = 3.09 ± 0.03. Most likely, 
the space density of optically-thick Lyman limit systems 
with N{HI) > lO^^cm"^ is enhanced within the red- 
shift range of the protocluster, tracking the galaxy den- 
sity. Our simulations indicate that the attenuation in the 
880-910 A range is more than twice as high for sightlines 
with intervening Lyman limit systems within one Lyman- 
continumn mean free path of the emission rcdshift, as 



opposed to sightlines without such absorbers. Therefore, 
we expect a corresponding enhancement in the average 
rest-frame 880 — 910 A attenuation for a composite spec- 
trum dominated by protocluster galaxies. The magni- 
tude of the enhancement in the average attenuation is 
not obvious, however, as demonstrated by the cases of 
our two individual Lyman-continuum detections. The 
sightline of C49 {z = 3.15) passes through the overden- 
sity, for example, yet it features very little if any Lyman- 
continuum absorption at 3.06 < z < 3.12. The redshift 
of D3 {z = 3.07) places it within the protocluster, yet its 
spectrum also appears fairly transparent in the Lyman- 
continuum region. A more careful treatment account- 
ing for the relative spatial distributions of absorbers and 
galaxies is required to quantify the effect of the protoclus- 
ter on the average IGM opacity in the Lyman-continuum 
region for our sample, yet we can state qualitatively that 
the z = 3.06 average Lyman-continuum attenuation fac- 
tor of 2.6 likely represents an underestimate. 

The second, related, caveat stems from the assumption 
in the IGM simulations that each randomly-generated 
sightline is probing a volume that is statistically- 
independent from the volumes probed by all other ran- 
dom sightlines. Since our sample is contained within an 
area of 9 x 13 comoving Mpc^ on the sky, ^ 50% of the 
resulting set of galaxy sightline pairs are separated by 
less than 5.7 como ving Mpc, the auto-corr elation length 
of z ~ 3 LBGs (|Adelberger et al.ll2005|) . This auto- 
correlation length is relevant since it has been demon- 
strated that, at least at N{HI) > 1 x 10-^'^cm^, in- 
tervening absorbers are associated with galaxies and 
therefore have similar c lustering properties (^eidcl 199^ 
lAdelberger et alJl2003D . The IGM opacity along neigh- 
boring sightlines in our sample must then be correlated, 
so it is not strictly correct to treat the galaxy sightlines 
as statistically independent of each other during our sim- 
ulations. The above caveats should be borne in mind 
throughout the remainder of the discussion. 

4.2. Variation in The Emergent UV Spectrum 

It is very striking that C49 and D3 appear to have 
such significant Lyman-continuum fiuxes and low ratios 
of UV to Lyman-continuum flux densities, while the re- 
mainder of objects in the sample are characterized by 
/goo non-detections. If there is significant variance in the 
emergent far-UV spectral shape in LBGs, understand- 
ing the origin of the variance will be a crucial compo- 
nent of characterizing the reionization of the universe. 
Therefore, the apparent variance in /1500//900 must be 
quantified more concretely, and related to the range of 
emergent far-UV spectral shapes, corrected for IGM ab- 
sorption (i.e. /1500//900 

First, as the sample spans a range in /1500, it is im- 
portant to confirm the fraction of our sample for which 
the observed /1500//900 ratios of C49 and D3 are de- 
tectable at the 3cr level. Since C49, for example, has a 
1500 A flux density that is above the median in the sam- 
ple, it is possible to detect larger /1500//900 ratios in this 
object than in an object such as MD14, which has one 
of the faintest 1500 A flux densities in the sample. By 
comparing the 1500 A fluxes of all the objects, with the 
3cr pixel to pixel noise limits in the 880 — 910 A wave- 
length range, we find that the observed /1500//900 = 7.5 
for D3 could have been detected in all of the objects. 
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The UV-to-Lyman-continuum flux-density ratio in C49, 
/1500//900 — 12.7, could have been detected in all but 
the two faintest objects, D7 and MD14. In fact, we find 
that the median detectable flux-density ratio for individ- 
ual objects in the sample was /1500//900 — 20.6. The 
most sensitive limit was for C32, in which we could have 
detected /1500//900 = 65.9, while the least sensitive was 
for D7, in which we could only have detected a ratio of 
/1500//900 = 8.3. ^ We therefore conclude that the finite 
depth of our observations did not prevent us from detect- 
ing objects similar to D3 and C49. If additional objects 
with similar /1500//900 had been in the sample, we would 
have detected Lyman-continuum flux from them. 

Second, in order to determine the range of emergent 
UV-to-Lyman-continuum ratios in the sample, we must 
take into account the uncertainty in IGM opacity along 
the line of sight in the Lyman-continuum region. For 
the sake of simplicity, though, we first consider the case 
in which all objects in the sample are characterized by 
the same IGM correction factor, and therefore, that the 
relative values (or limits) in /1500//900 translate directly 
into the relative values on /1500//900 icorr- In this case, 
we find significant differences among the emergent fiux- 
density ratios. Six objects have emergent UV-to-Lyman- 
continuum ratios greater than twice the ratio in C49, 
while 11 objects have emergent flux-density ratios greater 
than twice that in D3. 

Of course, variance in /1500//900 does not necessarily 
imply the same level of variance in /i5oo//900,corr - In- 
qualitative terms, if C49 and D3 happen to have sig- 
nificantly more transparent lines of sight than the av- 
erage at z 3, while other objects in the sample with 
900 A non-detections have sightlines with significantly 
more than average attenuation in the Lyman contin- 
uum, then the non-detections could easily have equiv- 
alent /1500//900, corr values. At z 3, the dispersion in 
IGM opacity for individual sightlines is fairly large. Once 
the uncertainty in IGM opacity for individual sightlines 
is included, the constraints are considerably weakened 
on the relative emergent flux-density ratios of C49, D3, 
and individual Lyman-continuum non-detections in the 
sample. In spite of this fact, we can still place some 
constraints on the sample variance in /1500//900 jcom 
by comparing C49 and D3 with the average of the re- 
maining 12 objects in the sample, for which the ob- 
served lower limit on the UV-to-Lyman-continuum ra- 
tio is /1500//900 > 43.0. This lower limit was estimated 
from the average spectrum of the objects without signif- 
icant Lyman-continuum detections, using the flux den- 
sity at 1500 A and the 3cr upper limit at 900 A, and 
includes both statistical and systematic uncertainties. ^ 
The uncertainty on the average IGM opacity for 12 ran- 
dom sightlines is signiflcantly smaller than for a single 
line of sight. For this 12-spectrum average, the 95.4% 

^ The ability to detect larger ratios of /1500//900 for D7 is lim- 
ited not only by the object's fainter-than-average flux at 1500 A, 
but also by its relatively low redshift (2 = 2.76), which causes 
rest-frame 880 — 910 A to fall at a bluer-than-average observed 
wavelength range with higher pixel-to-pixel noise than the typical 
value in our sample. 

The /1500//900 lower limit may be smal ler i f the systematic 
flux offset applied to the spectra (see section [3. II underestimates 
the true required value. However, we proceed with the assumption 
that the offset applied was the correct one. 
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Fig. 10 . — A zoomed-in comparison with the results of 
ISteidel et al.l 1200 1]V The composite spectrum from the current 
sample is plotted in black, while the Steidel ct al. 12001) composite 
is shown in r ed. The Lyma n limit is indicated with a vertical dot- 
ted line. The ISteidel et aTl (2)01) spectrum has been scaled to the 
level of the current composite at 1500 A. The pixcl-to-pixel noise 
level in the Lyman-continuum region of the Steidel ct al. (2005) 
curve is more than two times larger than in the new composite 
spectr um. Over t he region 880 — 910 A, the average flux density in 
the St eidel et al.l j200L,'l spectrum is greater than that in the new 
composite at high formal significance, despite the level of pixel-to- 
pixel noise. 



confidence interval on the correction factor for IGM at- 
tenuation yields /1500//900, corr > 13 — 20 (the factor 
13 [20] corresponds to the most [least] IGM attenua- 
tion) , though these numbers do not include the effects of 
the SSA22a galaxy overdensity on the Lyman-continuum 
opacity. As discussed above, the 95.4% confidence lim- 
its for C49 are /isooZ/goOjcorr = 0.3 — 7.5, while they are 
/i500 //900 icorr = 0.2 — 4.6 for D3. The lower limits on 
/1500//900 jcorr are unphysical, but here we are only inter- 
ested in the upper limits. Unless neglecting the SSA22a 
galaxy overdensity causes us to underestimate signifi- 
cantly the average IGM opacity for the 12 non-detections, 
we have indeed demonstrated that the emergent UV-to- 
Lyman-continuum ratios in C49 and D3 are significantly 
different from the average of the remaining objects in the 
sample with only Lyman-continuum non-detections. The 
challenge now is to explain this variation. 

5. IMPLICATIONS 

In the previous section, we found a sample average UV- 
to-Lyman-continuum fiux-density ratio, corrected for 
IGM absorption, of /1500//900 icorr = 22, which trans- 
lates to a relative Lyman-continuum escape fraction of 
fpscre i = 14%. The corresponding values in Steidel ct a^ 
l|2001|) are /i50o//900,corr = 4.8 (slightly different from 
the value of 4.6 presented in that paper, because of the 
different method used for correcting for IGM absorption), 
and fesc,rei — 63%. Therefore, using the average IGM 
opacities for the average emission redshifts for each sam- 
ple, we find that the estimate of fesc,rei for the current 
L BG sample is ~ 4.5 t imes lower than that of the sample 
in ISteidel et"an l)200H) . If we consider the 95.4% confi- 
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dence intervals on IGM opacity for a sample of 29 sight- 
lines at (z) — 3.40 and 14 sightlines at (z) = 3.06, we 
still find that fesc,rei for the present sample is at least 
3 times lower than that for the fSteidel ct al. (200]J) 
sample. Figure lTHl shows a zoomed-in view of the Lyman- 
continuum regions for the .Steidel et al.i (;2001,'l composite 
(in red) and the new composite (in black), with the spec- 
tra scaled to the same value at 1500 A. The'Steide l et al.1 
l|2"001,) spectrum, while noisier, is still characterized by 
a mean flux-density value in the 880 — 910 A region 
that is more than three standard deviations (as measured 
by pixel-to-pixel fluctuations) higher than the mean flux 
density in our new LBG composite spectrum. 

If we adopt the new estimates for (/1500//900 jcorr ) and 
{fesc,rei) SlS representative of the LBG populati o n, it i s 
possible to scale the values from iSteidel et al.1 l)200lD . 
to estimate the global contributions of LBGs to the 
ionizing background at z ~ 3. For this estimate, 
we use (/i5oo //900 icorr ) to convert the published LBG 
1500 A luminosity function l|Adelberger fc Steideill200f]() 
to one at 900 A. This 900 A luminosity function is in- 
tegrated down to O.IL* to estimate the comoving emis- 
sivity from LBGs at the Lyman limit. Taking into ac- 
count the mean free path of Lyman-continuum photons 
at z - 3, which is Az ~ 0.18, we find Ji.(900) - 
2.6 X 10~^^ erg s~^cm~ ^Hz~^sr~^. ^ If the composite 
spectrum in ISteidel et~a l. ( 2001) is adopted as charac- 
teristic of the entire LBG population, the inferred con- 
tribution to the ionizing background is correspondingly 
higher: J,.(900) - 1.2 x 10"^^ erg s~^cm-^Hz~- ^sr-^ 
Because the galaxies included in the ISteidel et alJ com- 
posite spectrum are drawn only from the bluest quartile 
of LBGs, that spectrum may only be representative of 
~ 25% of LBGs. Under the assumption that the bluest 
quartile of LBGs are also the most likely to have es- 
caping Lyman-continuum emission and that the remain- 
ing, redder, LBG population is charact erized by a neg- 
ligible escape fraction, the Stei del et al.l composite spec- 
trum translates into a Jy(900) value consistent with the 
one presented here. Combining a new measurement of 
the faint end of the z ~ 3 QSO UV luminosity function 
with existing data on the bright end, and assuming 100% 
escape fraction of ionizin g radiation at all QSO lumi- 
nosities, l^n^OD present estimates of the QSO 
contribution to the ionizing background at z ~ 3. This 
estimate is Ji,(900) ~ 2.4 x 10"^^ erg s~^cm~^Hz"^sr"\ 
very similar to the estimate for LBGs presented in this 
work. ^ If galaxies and QSOs contribute roughly equal 
amounts to the ionizing background at z ~ 3, this should 
be reflected in the spectrum of the metagalactic ion- 
izing radiation field, which in turn determines the rel- 

^ The Steidel et al. (2001) result was derived using an Einstein- 
de Sitter cosmology. However, Jj,, which is directly proportional to 
the product of comoving emissivity and Lyman-continuum mean 
free path, is roughly independent of cosmol ogy. Therefor e, we are 
justified in applying a simple scaling to the ISteidel et all value for 
J;. (900). 

® The Hunt et alj 120041) work also assumed an Einstein-de Sit- 
ter cosm ology, and therefore makes for a fair comparison with the 
ISteidel et al.i Ii200ll) results, and, by extension, the current estimate 
of Ju from galaxies. Because the estimate of Ji, is roughly inde- 
pendent of cosmology, the values presented for QSOs and galaxies 
in Hunt et al. 1 2004) and the current work, respectively, can also 
bo compared consistently with Ju values ba sed on an f2m = 0.3, 
= 0.7 cosmology (e.g., those in lBolton et al. 2005). 



ative ra tes of ionization of inte r galactic hydrogen and 
helium llHaardt fc MadaiH 119961: ICen fc Ostrikea 120061: 
iBolton et al.ll2006D . 

It is valuable to compare our results with recent esti- 
mates of the metagalactic hydrogen ionization rate de- 
rived using independent methods. Modeling the ob- 
served Lya forest opacity in QSO absorption spectra 
using a suite of high-resolution hydrodynamical simula- 
tions, ^olto^^ al. (2005i) conclude that the global hy- 
drogen ionization rate at z ~ 3 is Fjjj = 8.6 ± 3.0 x 
lO^^'^s^-'^. W ith a spectral slope for J(z/) of ly^^-^, as in 
iMadau et all ((1999), this Fjjj corresponds to Ji,(900) 
3.4 X 10~^^ erg s~-^cm~^Hz~^sr~^. If the shape of the 
ionizing spectrum lies between and v'^, the con- 
version between Fjjj and J^{900) ranges between 4.2 
and 2.2 x 10~^^ erg s~^cm~^IIz^^sr~^. Adopting the 
ISteidel et~al . (2001) composite spectrum as representa- 
tive of the entire LBG population therefore translates 
into a LBG contribution to the ionizing background ex- 
ceeding by more than a factor of ~ 3 the value of 
Jiy(900) inferred from modeling the Lya forest absorp- 
tion. In contrast, our new estimate for the contribu- 
tion from LBGs to the ionizing background does not ex- 
ceed the Lya forest-inferred result, and, furthermore, the 
new sum of the c ontributions from galaxies and QSOs 
fHunt et al.'2004') is roughly consistent with the total H I 
ionizing radiation field presented in Bolton ct al. (2005). 
It is worth not ing that the quasar proximity effect study 
of IScott et alJ (^2000.) indicates a value for the ionizing 
backgroun d at z ^ 3 th at is roughly twice as large as 
the one in IBolton et al.l and is arguably a more direct 
method. The uncertainty in the quasar proximity effect 
result is large enough that the two methods still yield 
values for Jy(900) that are actually statistically consis- 
tent with each other, but it will be important to reconcile 
the systematic offset in the mean inferred J,y(900) from 
these different techniques. 

Estimating the global contribution from LBGs to the 
ionizing background is based on the assumption that we 
have determined the average shape of LBG spectra be- 
tween the UV and Lyman-continuum regions. This as- 
sumption should be viewed with a degree of skepticism 
for more than one reason. First, for the two samples of 
LBGs with flux measurements in the Lyman-continuum 
region, two different values of /1500//900 icorr have been 
derived. This fact is not particularly surprising, given 
that the two samples are characterized by different aver- 
age rest-frame UV spectroscopic properties, which may 
in t urn be related to the emergent far-UV spectral shape. 
The lSteidel et al.l lj200m sample is characterized on aver- 
age by stronger Lya emission, weaker low-ionization in- 
terstellar absorption lines, and less dust extinction. The 
weaker interstellar absorption and dust extinction may 
both be correlated with hig her escape fra c tion o f Lyman- 
continuum photons in the ISteidel et al.l l)2001(l sample. 
Furthermore, neither of the two Lyman-continuum sam- 
ples is completely representative of the total LB G spec- 
troscopic sample. Both the current and iSteidel et alJ 
samples have brighter than average UV luminosities - 
most striking in the current sample, in which all the 
galaxies are brighter than L* of the LBG UV luminosity 
function. While the current sample has an average in- 
ferred dust extinction very close to that of the total LBG 
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sample, the average Lya emission and interstellar absorp- 
tion line strengths are, respectively, weaker and stronger 
than in the total LBG sample. While the lSteidel et alj 
sample has similar average Lya emission and interstellar 
absorption properties to the total LBG sample, it is sig- 
nificantly bluer in the UV continuum. Furthermore, most 
of the galaxies in the current sample reside in a signifi- 
cant overdensity, which represents an atypical large-scale 
environment and may affect the level of IGM opacity in 
the Lyman-continuum region. Therefore, we may have 
only succeeded so far in measuring /1500//900 jcorr for dif- 
ferent types of LBGs drawn from the bright end of the 
luminosity function, some of which reside in a special en- 
vironment, without much constraint on the average ratio 
for the total LBG sample. 

The second reason for skepticism is that, even within 
our sample, we have demonstrated a significant varia- 
tion in fi5oo/ f 900, corr- C49 and D3, the two galaxies 
for which we detected flux in the Lyman-continuum re- 
gion, must have significantly different /1500//900, corr ra- 
tios from at least the average of the remaining galaxies 
in our sample. The rest-frame UV properties of C49 
and D3 do not immediately stand out, however, relative 
to those of the objects without Lyman-continuum detec- 
tions. Furthermore, these two galaxies have properties 
strikingly different from each other! Since we do not 
yet have a coherent physical explanation - in terms of 
stellar population age or ISM geometry - for why these 
two galaxies alone exhibit significant flux in the Lyman- 
continuum region, we are unable to show that our small 
sample is representative even of bright LBGs. As the 
intrinsic {L i^oo /L goo) int ratio depends so sensitively on 
stellar population age, it is important to obtain near- and 
mid-IR data, which will allow us to constrain the stellar 
populations of these galaxies. On the other hand, the 
variation in /1500//900 icorr among our sample may sim- 
ply reflect the fact that the escape of Lyman-continuum 
emission is an anisotropic process, resulting in a complex 
spatial distribution of Lyman-continuum emission rela- 
tive to that of the UV-continuum emission determining 
the positions of our spectroscopic slits. 

What is required is a statistical sample of LBGs with 
deep Lyman-continuum measurements. This sample 
should be roughly an order of magnitude larger than 



the current one, span a larger range in luminosity, with 
multi-wavelength data from which constraints on stel- 
lar population ages can be obtained. Furthermore, to 
enable robust correction for IGM opacity along the line 
of sight, the future sample must contain objects prob- 
ing at least several statistically independent cosmic vol- 
umes - as opposed to the objects in this work, all of 
which are contained within 9 x 13 comoving Mpc^ on the 
sky and most within a previ ously-discovered galaxy over- 
density ijSteidel et al.l[l99 8lh. This larger sample should 
also contain a more significant set of Lyman-continuum 
detections, whose properties can be characterized on 
a statistical basis, relative to those of objects without 
detections. A lack of intrinsic difference between the 
two populations would support the idea that variation 
in /1500//900, corr reflects anisotropic Lyman-continuum 
emission being sampled with different degrees of effi- 
ciency. Only with a survey satisfying the above require- 
ments will we be truly justified in translating the LBG 
UV luminosity function into a constraint on the global 
contribution of LBGs to the ionizing background. How- 
ever, what we have shown here, the first detection of ion- 
izing radiation from individual LBGs at 2; ~ 3, represents 
significant progress in the right direction. 
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TABLE 1 

LRIS Observations of Bright Lyman-Break Galaxies in the SSA22A Field 



Object 


R.A. (J2000) 


Dec. (J2000) 




G -TV^ 


Un - G^ 


ZerJ' 


2^ 


LRIS-B 

Exposure (s)*^ 


LRIS-R 

Exposure (s) 


D17 


22:17:18.88 


00:18:16.89 


24.27 


0.45 


2.01 


3.0851,3.0945 


3.0697 


29200/69040 


61000 


C24 


22:17:18.94 


00:14:45.65 


23.86 


0.78 


> 2.77 


3.1026 


3.0960 


29200/79240 


61000 


C49 


22:17:19.81 


00:18:18.64 


23.85 


0.59 


> 3.04 


3.1601,3.1656 


3.1492 


29200/58040 


61000 


C35 


22:17:20.23 


00:16:52.33 


24.18 


0.95 


> 2.52 




3.0980 


29200/79240 


61000 


C47 


22:17:20.23 


00:17:32.45 


23.84 


0.60 


> 3.14 


3.0748 


3.0651 


29200/79240 


61000 


C41 


22:17:24.45 


00:17:14.75 


23.80 


0.18 


> 3.49 


3.0236 


3.0166 


29200/79240 


61000 


C32 


22:17:25.62 


00:16:13.10 


23.68 


0.67 


> 3.19 


3.2895,3.2991 


3.2921 


29200/79240 


61000 


Cll 


22:17:25.68 


00:12:35.30 


24.20 


0.47 


> 2.95 


3.1014 


3.0962 


29200/79240 


61000 


MD23 


22:17:28.01 


00:14:29.73 


24.14 


0.46 


1.92 


3.0845 


3.0753 


29200/79240 


61000 


D7 


22:17:28.79 


00:12:47.07 


23.50 


0.62 


2.14 




2.7564 


29200/47840 


61000 


C16 


22:17:31.96 


00:13:16.11 


23.64 


0.98 


> 2.88 




3.0651 


29200/79240 


61000 


D3-apl 


22:17:32.42 


00:11:32.97 


23.37 


0.97 


2.58 


3.0716 


3.0649 


29200/47840 


29200 


D3-ap2 


22:17:32.42 


00:11:32.97 


23.37 


0.97 


2.58 


3.0767 


3.0677 


29200/47840 


29200 


MD14 


22:17:37.93 


00:13:44.21 


24.14 


0.86 


2.25 




3.0969 


29200/29200 


29200 


C22 


22:17:38.18 


00:14:32.06 


24.46 


0.54 


> 2.72 




2.8765 


29200/58040 


61000 



'^AB magnitude system; effective wavelengths are 4730, and 6830 A for G and Tl, respectively (see also lSteidel et al.ll200^ . 

''Spectroscopic redshift measured from the observed wavelength of Lya emission. Three objects, D17, C49, and C32, exhibit double-peaked Lya profiles. For these objects, we list the 
emission redshifts corresponding to both Lya peaks. 
'^Average spectroscopic redshift measured from the observed wavelengths of strong low-ionization interstellar absorption lines. 
LRIS-B exposure times are listed for A < 4000 A (left) and 4000 A< A < 5000 A (right). 
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TABLE 2 

Properties of Different LBG Samples 



This Work Total LBG Sample Steidel et al. (2001^ Sample 

N I 14 811 29 

(z) 3.06 ±0.09 2.96 ±0.29 3.40 ±0.18 

(n) 23.92 ±0.32 24.60 ±0.56 24.34 ± 0.38 

{Muv)"^ -21.63 ±0.32 -20.88 ±0.58 -21.40 ± 0.38 

{G-n)^ 0.65 ±0.23 0.61 ±0.29 0.89 ±0.18 

(EiB-V))" 0.11 ±0.08 0.13 ±0.09 0.07 ± 0.05 

(VKo{LyQ)><i 1.7 A 15.2 A 14.1 A 

(WollS))'^ -2.2 A -1.8 A -1.7 A 



^AB absolute UV magnitude, derived from the apparent H magnitude. At 2 ~ 3.0 — 3.4, the 7^-band effective wavelength corresponds to 
1550 - 1700 A. 

''The larger < G — R > of the lSteidel et all 120011) sample reflects the increased average IGM opacity at z = 3.40 relative to 2 = 3.06 and 
2 = 2.96. In fact, this sample appears to have bluer intrinsic colors on average than the two other samples. 

•^Dust extinction, parametrized by E{B—V), is derived from the observed G—R color corrected for the effects of intergal actic absorption and 
intrinsic Lyo emission/absorption. The corrected G — R color is modeled, and E{B — y) inferre d, using an unreddened teruzua^^^hariort 
jE003) 300 Myr contirmous star-formation model template spectrum and a lCalzetti et al.l 12003) attenuation curve (see lSteidel et alj|199a 
lAdelbergcr k. Steidel 2000; Shaplev et al. 2001). 

Rest-frame Lyo and low-ionization interstellar absorption equivalent widths have been measured from the average composite spectrum 
for each sample. Positive numbers refer to emission and negative numbers to absorption. In the case of Lya, t he reported equival ent width 
represents the total profile, including both emission and absorption components. For both the Total LBG and Steid el et alj 12001]) samples, 
the total and emission profiles have very similar equivalent widths. In the case of the current sample, the total equivalent width is close 
to ^ A while the emission equivalent width is ~ 5 A rest. The reported low-ionization interstellar absorption equivalent width for each 
sample represents the average of the four strongest features: Si II A1260, O I A1302 + Si II A1304, C II A1334, and Si II A1526. 



